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ABSTRACT 

Modern population synthesis models estimate that 50% of the restframe if -band light is produced 
by TP-AGB stars during the first Gyr of a stellar population, with a substantial fraction continuing to 
be produced by the TP-AGB over a Hubble time. Between 0.2 and 1.5 Gyr, intermediate mass stars 
evolve into TP-AGB C stars which, due to significant amounts of circumstellar dust, emit half their 
energy in the mid-IR. We combine these results using published mid-IR colors of Galactic TP-AGB 
M and C stars to construct simple models for exploring the contribution of the TP-AGB to 24/^m 
data as a function of stellar population age. We compare these empirical models with an ensemble 
of galaxies in the CDFS from z — to z — 2, and with high quality imaging in M81. Within the 
uncertainties, the TP-AGB appears responsible for a substantial fraction of the mid-IR luminosities 
of galaxies from z = to z = 2, the maximum redshift to which we can test our hypothesis, while, 
at the same time, our models reproduce much of the detailed structure observed in mid-IR imaging 
of M81. The mid-IR is a good diagnostic of star formation over timescales of ~ 1.5 Gyr, but this 
implies that on-going star formation rates at z = 1 may be overestimated by factors of ~ 1.5 — 6, 
depending on the nature of star formation events. Our results, if confirmed through subsequent work, 
have strong implications for the star formation rate density of the universe and the growth of stellar 
mass over time. 

Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: stellar content — infrared: 
galaxies 



1. INTRODUCTION 

The infrared provides a critical window into obscured 
star formation in our Galaxy and in nearby galaxies 
(e.g. Helou et al. 1988 and many since), allowing us to 



peer into a range o fstar forming environments at the 



prese nt epoch (e.g. Roussel et al. 2001 Calzetti et al 



2007 1 . Space-based IR observations have improved and 



extended our vi ews to greater dist ances and large look- 
back times (e.g. Salim et al. 2009), making these wave- 
lengths crucial for studies ot galaxy str ucture (see, e.g., 
Regan et al. |2004 ) and evolution (e.g. Papovich et al. 
2007p . I'he 24/um data, in particular, provided the first 
of what were to be unbiased histories of the star forma- 



tio n rate density of the universe over time (Le Floc'h et 



al. 



2005). 



Such analyses utilize template SE Ds in order to esti- 



mate IR bolometric corrections (e.g. 


Chary & Elbaz 2001 


Dale & Helou 


2002 Rieke et al. 


2009). These, in turn 



are used with empirical calibrations oi on-going SFRs de 
rived from nearby galaxies under the assumption that the 
mid-IR flux arises from the reprocessed light from young, 
luminous stars. However, these circumstances are not 
well understood, partly because of the different physical 
mechanisms and timescales probed by the IR compared 
to recombination lines or UV emission (e.g. Kennicutt 



1998 Salim et al. 



2009) 



Here we recognize that when intermediate mass stars 
join the Thermally-Pulsating Asymptotic Giant Branch 
(TP-AGB), they do not uniformly have the colors of 
lower-mass TP-AGB M stars but, as C stars, are partic- 
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ularly luminous in the mid-IR, with m ost of their energy 
emitt ed between 20fim and 45/^m (e.g. Guandalini et al. 
2006|. Knowing the imp ortance of the TP-AGB at rec 



and near-IR wavele ngths ( Maraston|2005| [Bruzual| 2009 



Conroy et al.||2009 ) , and, using near- to mid-IR colors o 
Galactic TP-AGB populations, we empirically calibrate 
the contributions of such stars to the integrated mid-IR 
luminosities of stellar populations. Using simple mod- 
els, we show that the mid-IR luminosities of galaxies are 
specifically in the sensitivity of the mid-IR to the amount 
of stellar mass formed in the previous 1.5 Gyr, naturally 
complimenting the optical and near-UV. 

2. OBSERVED PROPERTIES OF THE TP-AGB 

Our calculations explicitly rely on the fraction of K- 
band light produced by the TP-AGB as a function of the 
age of a simple stellar populati on (SSP). Fig ure 13 of 
Maraston| ( f2005| ) and Figure 1 of |Bruzual| j2009[ ) show the 
TP-AGB s contribution to the A -band as a timction of 
time, here approximated in our Figure [Ha). Between 0.1 
and 1.5 Gyrs, the TP-AGB is responsible for roughly half 
the luminosity at K . At solar metallicity, and between 
0.2 and 1.5 Gyrs, stars have sufficient mass to becom e 



C stars (Marigo & Girardi 

At earlier times hot bott om burning prevents sta rs trom 



2007| |Marigo et aT^|2008| 



inverting C/O ratios (e.g. Marigo & Girardi 2007[ ). Lower 
mass stars that contribute to the A -band over a Hubble 
time also remain O-rich. Consequently, mo st assume the 
TP-A GB to have the color of M stars (e.g. |Willner et~aL 
2004| . 

Knowing how the luminosity of the TP-AGB evolves at 
K allows one to compute the integrated mid-IR luminos- 
ity if the near- to mid-IR color for the ensemble of stars 
on the TP-AGB is known as a function of time. There 
now exist sufficient data in the literature for Galactic 
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Figure 1. The mid-IR contributions of the TP-AGB. (a) Simplified schematic for the fraction of iC-band light produced by the TP-AGB 
for an SSP. This approximation allows us, in tandem with observed near- to mid-IR colors of TP-AGB stars, allows us to estimate mid-IR 
luminosities. The blue region specifics the time over which stars can become C rich and emit half their energy in the mid-IR. Outside of 
this age range, we use the ensemble color of the Galactic TP-AGB M stars, (b) The correlation between £24/(1+2) an d Lk f° r galaxies in 
CDFS at 2 ~ 1 (courtesy of |Wuyts et al,|2008| l, where those with 24/^m flux densities significant at a level less than 2a are shown in black, 
less than 5<r in green, and the rest m orange. Non-detections are shown as red points along the bottom of the figure. The maximum and 
minimum mid-IR, luminosities of galaxies are defined by the ensemble colors of TP-AGB stars, shown by the blue and red lines, for young 
and old stel lar populations, respectively , (c) Galaxies in CDFS at z ~ 1.8. Lavender circles mark passively evolving galaxies at z ~ 2 in 
the HUDF < |Maraston et a.l.||2006] |2007[ |. (d) Galaxies in CDFS and the HUDF where the blue and red circles denote the colors of the 
TP-AGB C and 1V1 ensembles in each of the four bands. The solid line traces the MIPS 24/xm bandpass where fc-corrections have been 
computed using the |Dale fc Helouj p002| templates. 



TP-AGB stars to constrain such colors, and we adopt 



the mid-IR bands 



work o f Guandalini et al. (|2006|) and Guandalini & Busso 
(120081). 



14.61 . and [21.3] from the 



Using 200 Gal actic TP-AGB C st ars from |Guandahm| 
et al. ( |2006[ ) and Busso et al. (2007), and conservatively 
excluding stars redder than K — [8.8] > 10 mag to re- 
duce any bias towards IR-luminous post-AGB objects, 
we compute mean colors for the population of dust- 
enshrouded Galactic TP-AGB C stars of (K - X) c = 
5.66 ±1.16, 6.41 ±1.15, 6.31 ± 1.06, and 6.87 ± 0.86 mag 
for the four ba ndpasses, respectively ( Vega magnitudes; 
zeropoints from |Guandalini et al.|2006| . We do not claim 
that all stars with sufficient mass have these colors upon 
joining the TP-AGB, but that they are dusty and IR lu- 
minous, such that when averaged over ti mescales longer 



than the T P-AGB lifetimes (10 6 yrs; e.g. |Marigo fc Gi 
rardi 2007) this color is a good approximati on tor the 
ensemble. Using Guandalini & Busso (|2008 ), we com- 
pute ensemble colors of {K - X) M = 1.17±~0.24, 2.32 ± 
0.34, 1.61 ± 0.41, and 2.88 ± 0.58 mag for TP-AGB M 
stars. 

Formal errors in these colors are smaller than the 



systematic uncertainties, which have been listed above. 
These were computed by varying the color bounds of 
the TP-AGB and by using several different estimators. 
These uncertainties dominate the errors in subsequent 
modeling. TP-AGB ensemble colors are critical, and use 
of the Galactic samples assumes that these stars ade- 
quately sample the phases of TP-AGB evolution. 

We now combine Figure Ilia) with the knowledge that 
the C star ensemble color should be used for SSP ages 
within 0.2 < t < 1.5 Gyr, and the M ensemble color 
for all other ages. Whether stars have sufficient mass to 
become C stars depends on mass-loss and the efficiency 
of dredge-up, both of wh ich are sensitive to metallicity 
(Marigo & Girardi 2007). For the purposes of making 
initial estimates of the contribution of the TP-AGB in 
the mid-IR, we ignore these important effects. 

3. IMPLICATIONS FOR GALAXIES 

The colors given in the previous section imply that for 
ages less than 1.5 Gyr, the TP-AGB emits more than 
twice as much energy in the [12.5] band as it does in 
K. In other words, a young galaxy at z = 1 with a 
Lk = lO 11 !^,© should have a luminosity of £i2.5^ m = 
1.6 x IO^LqHz" 1 , equivalent to a LIRG with L T ir — 
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Figur e 2. (a-c) Bicolor diagrams u z — g z vs g z — z z , K z — [24] vs u z — g z , and u z — [24] vs g z — z z for galaxies at 0.7 < z < 1.1. Galaxies 
from |Wuyts et al.| ( |2008[ > are also shown, cut at L > 10 Lq in the restframe V-band for clarity. Blue circles show the simple models when 
logSSlK > — 10, and red triangles when logSSFR < —10. The red vector illustrates the effect of Ay = 1 mag of extinction of starlight 
< |Calzetti et al.|2000[ l, with no reprocessing of UV light into the mid-IR. The u z g z z z K z denote redshifted bandpasscs. In (d) we plot the 
correlation between model 24/xm luminosities to star formation rates computed within the 1.5 Gyr window o ver which TP-AGB C stars 
are visible. The violet dashed line is the correlation of 12/xm luminosity with SFR from|Chary Sz Elbaz ( 2001[ |. 

1.3 x 10 n Lp, ( [Chary fc Elbai]|2001| ). 

Figures [ijb) and (c) illustrate this key point at red- 
shifts where the MIPS 24/im band corresponds to [12.5] 
and [8.81. When stellar populations are producing TP- 



AGB C stars and their fraction of if-band light is 50%, 
one expects the maximum output in the mid-IR relative 
to K, shown by the blue line. Very young stars (< 200 
Myr) and those older than 1.5 Gyr will produce only 
TP-AGB M stars, delineating the minimum output in 
the mid-IR relative to K, shown by the red li ne. These 



predictio ns are compared t o galaxies in CDFS dWuyts et 
al] [2008l ) and the HUDF ( |Maraston et^[2^0rJt - pD07| ; 



where restframe magnitudes were derived by interpolat- 
ing the published photometry. In Figure |T|(cl) we expand 
these simple comparisons to the redshift range z = to 
z = 2. 

It is striking that the circumstellar material of evolved 
stars alone, up to 1.5 Gyr old, can explain the observed 
mid-IR luminosities of galaxies back to z = 2. For galax- 
ies near the lower bound, a substantial fraction of their 
24/im fluxes may originate purely from old stars, even 
at early times, where the actual fractions will depend on 
more precise modeli ng. This conclusion is fu lly consistent 
with earlier work by Maraston et al. (20061, whose mod- 
eling of passively evolving galaxies in the H LJDF (see Fig. 
[l]) at z ~ 2 highlighted contributions of the TP-AGB 
at restframe near-IR wavelengths. Using the [8.8]-[21.3] 



colors of the TP-AGB ensembles derived above, these 
galaxies should have modest detections of 1 — 10/xJy in 
future observations at 60/im. 

In the next section we construct simple star-formation 
histories to study the observed correlations of the mid-IR 
with restframe optical bands. By doing so we can check 
our hypothesis against the well-observed evolution of the 
global properties of galaxies. More sophisticated efforts 
to incorporate the TP-AGB int o population syn thesis 
have a number of unsolved issues ( Conroy et al. 1 20091 ) and 
physically motivated models Marigo et al.| ( |20U8| F"cTo not 
yet adequately match the broa d range ot observed mid - 
IR colors of TP-AGB C stars (|Guandalini et al ||2006| ). 
Therefore we explore the mid-IR luminosities of stellar 
populations as a function of redshift using our empirical 
calibration of the mid-IR luminosity of the TP-AGB. 

4. SIMPLE FORMATION HISTORIES AND THE MID-IR 

The ensemble colors of TP-AGB stars have been shown 
to describe the range of mid-IR luminosities of normal 
galaxies back to z = 2. We now use model star forma- 
tion histories (SFHs) to compare the evolution of stellar 
populations in the mid-IR w i th the evolution o f galaxy 
colors . We use the Maraston (20051 models, the |Kroupa 
( [2001] IMF, solar metallicityT^TP-AGB X-band iu- 
mmosity fraction that evolves as in Figure [T] and the 
mid-IR C and M ensemble colors given above. No dust 
components outside of those implicit in the colors of the 
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Table 1 

Star Formation History Parameters 



Type 


Mstars,z=0 


20 




T 


N T 


Red 


3 x 10 11 


5 


1.0 


1 x 10 9 


3 


Red 


8 x 10 10 


5 


2.0 


5 x 10 8 


2 


Red 


6 x 10 10 


5 


2.0 


1 x 10 9 


oo 


Red 


4 x 10 10 


5 


1.5 


1 x 10 9 


oo 


Red 


1 x 10 11 


5 


1.0 


5 x 10 8 


3 


Red 


9.5 x 10 10 


5 


4.5 


3 x 10 8 


2 


Blue 


3 x 10 10 


5 


1.8 


3 x 10 9 


oo 


Blue 


5 x 10 10 


5 


1.5 


2 x 10 9 


oo 


Blue 


2 x 10 11 


5 


2.5 


1.5 x 10 9 


oo 


Blue 


1.5 x 10 11 


5 


2.5 


2.0 x 10 9 


oo 


Blue 


2 x 10 10 


5 


1.2 


5.0 x 10 9 


oo 


Blue 


2.5 x 10 10 


5 


0.8 


5.0 x 10 9 


oo 



Galactic TP-AGB stars have been added. 

Several SFHs have been constructed, with parame- 
ters given in Table [I] We use exponentially declining 
SFHs, or T-models, modified by having each model be- 
gin at z Q = 5, with an exponential rise of 1 Gyr until a 
peak at redshift z p . The models explicitly terminate N t 
timescales after z p . The models have been normalized to 
stellar masses ranging from 2 x 1O 1O M and 3 x lO n M0 
at z = 0. Parameters were chosen to produce colors con- 
sistent with a diversity of blue and red galaxies at the 
present epoch. 

4.1. Galaxy Colors at z ~ 1 

In Figures [^a-c) we compare the colors of the mod- 
els to the galaxies in CDFS between 0.7 < z < 1.1 
(Wuyts et al. 20081 using the bicolor diagrams u z — g z 
vs g z - z z , K z - [24] vs u z - g Zl and u z - [24] vs g z ~ z z , 
(u z g z z z K z denote redshifted passbands). Blue triangles 
indicate times when logSSFR > —10, and red trian- 
gles when logSSFR < —10, where SFR refers to on- 
going star fo rmation rates. The extinction of starlight 



by Ay = 1 (Calzetti et al. 2000[ ) is shown by the red 
arrow. With no additional tuning, these SFHs mimic 
both the blue star forming sequence and the red qui- 
escent sequence, while simultaneously reproducing the 
correlations of 24/im emission with optical colors. 

Figure ^td) shows that the simple models reproduce 
the correlation b etwee n 12/im luminosity and SFRs of 
Chary fc Elbazj (20011, but only when we average the 
SI R over the 1.5 Gyr prior to the epoch(s) of observa- 



tion. |Chary fc Elbaz (2001 ), and others, have calibrated 
the underlying SFK as a function of mid-IR luminosity 
and, though our models reproduce this correlation, the 
timescale probed by this relation is much longer than has 
been assumed in the past. In other words, the mid-IR is 
explicitly sensitive to star formation over the timescales 
that galaxies produce populations of TP-A GB C stars, 
or 1.5 Gyr (confirming the earlier results of Salim et al. 
2009| . 

Together these results also imply that the relation- 
ship between rates of on-going star formation and mid-IR 
fluxes will be complicated by the detailed SFHs within 
a 1.5 Gyr windows. Bursts that occur at z = 1.5 and 
decay with r = 1 Gyr, still produce substantial mid-IR 
fluxes at z = 1, resulting in estimated SFRs 2.3x higher 
than actual on-going rates. For r = 0.5 Gyr, the factor 
is 6.4. Using the ~ 3 Gyr exponential timescale for the 



decline of the star formation rate density of the universe 
at late times implies that the mid-IR is an over-estimate 
by at least 50%. Taken a step further, post-starburst 
and "green valley" galaxies may simply be luminous in 
the mid-IR because TP-AGB C stars continue to appear 
long after the cessation of major star formation activity. 

Calculating the impact on the star formation rate den- 
sity ultimately requires knowing the frequency, duration, 
and intensity of major mass building events. Given the 
unknown duty cycles of major star forming events, mid- 
IR-based star formation rates of individual galaxies at 
z = 1 should be treated with skepticism at a level of at 
least a factor of two, with additional uncertainties in star 
formation rate densities. 

4.2. Normal galaxies at Late Times: M81 

Our models have strong implications for the origins 
of mid-IR emission in nearby, resolved systems. Here we 
compare data on the nearby galaxy M81 with our model's 
implied correlations of the mid-IR with (1) luminosity at 
K, and (2) the relative amount of star formation in the 
previous 1.5 Gyr of evolution. 

SDSS u and 2MASS K s images are shown in Figures 
[3]ja) and (b). Binning these to increase the S/N ratios 
per data point, we show color-color diagrams for points 
in the galaxy down to [Ik = 20 mag per arcsec 2 in (c) and 
(d). Our smooth SFHs at late times are shown by the 
thick black lines. Colors indicate galaxy radius: R < 1.5 
min (red), 1.5 < R < 3 min (green), 3 < R < 5 min 
(blue), and R > 5 min (violet). For simplicity, we assume 
points only fall redward of these models in u — g vs g — z 
because of extinction and thus derive a map of Ay across 
M81. These extinctions allow us to correct u—g and infer 
the ratio S24/ 1 Sk at each position in the galaxy using Fig. 
Ic). 

We show the observed 24/^m map along with the re- 
sulting model in Figures Qa) and [4jb) . Aside from the 
limiting S/N of the u and K s data, the agreement is 
good. We also find similar results in 8/zm but have not 
included them owing to space considerations. In Fig. 
4](c) we show the encircled flux densities in both 24/^m 
and 8/im, and one can directly see that the luminosity 
and structure of M81 is consistent with TP-AGB ori- 
gins. Within R < 3 arcmin, the simple model is missing 
~ 20% of the flux, but given uncertainties in using u — g 
color as a proxy for age, uncertain extinction corrections, 
and unaccounted metallicity effects, we cannot rule out 
additional contributions from interstellar cirrus. In re- 
gions of intense star formation, which are known to be 
dust enshrouded, the model also underproduces the mid- 
IR fluxes, but such regions do not comprise the bulk of 
M81's luminosity. 

We conclude that the extended, diffu se mid-IR emis- 
sion noted by Kennicutt et al. (20091 arises from the 
stellar populations of the galaxy. This result is consis - 







tent with the analysis of M33 by Verley et a"T] ( 2009 1 
who employed a more limited set of models and argued 
that dusty circumstellar envelopes of evolved AGB stars 
may be the source of the bulk of M33's diffuse 8/xm and 
24/xm emission. 

5. SUMMARY 

The contribution of the TP-AGB to the if-band has 
been combined with with the mean IR colors of Galactic 
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Figure 3. M81 in the optical and near-IR. (a) SDSS M-band image, (b) 2MASS i^s-band image. In both cases we show circular radii of 
1.5,3.0,5.0, and 7.0 min. (c) The ratio of 24^tm to K a flux density vs u — g color. The colors of the points indicate spans of radius with 
red (R < 1.5 min), green (1.5 < R < 3 min), blue (3 < R < 5 min), and violet(i? > 5 min). The thick black line shows the locus of our 
smooth SFHs at late times, assuming no extinction of star light, (d) The u — g vs g — z bicolor di agram, with the sam e models shown by 
the thick black line. The red vectors illustrate the effect of Ay = 1 mag of extinction of starlight (Calzetti et al. 2000} . 

TP-AGB C and M stars in order to estimate the con- 
tributions of both young and old stellar populations to 
mid-IR observations of galaxies. Without tuning, we find 
that the resulting mid-IR luminosities of the TP-AGB 
can reproduce the MIPS 24/im fluxes for galaxies back 
to at least z — 2 in a manner consistent with restframe 
optical colors. We have also tested the validity of the 
model on local scales in the galaxy M81 and find rea- 
sonable agreement. The origins of correlations between 
optical colors and mid-IR luminosities seen by others, 



There is little doubt that star formation and the 
growth of stellar mass was occurring more rapidly in 
the distant universe than today, but the nature of that 
grow th has remained large ly unknown. Earlier results 



(e.g. Le Floc'h et al. 



|2005[ ) had implied that - 1/3 of 
the stellar mass at the present epoch was formed after 
z = 1 — a result that appears to be at odds wit h the 
evolution in the stellar mass function to z 



1 (e.g. Cira- 



such as Salim et al. ( 2009) , can now be understood. 

With careful modeling of SEDs from the UV through 
the mid-IR, more detailed histories of star formation 
should be possible. Unfortunately, stellar spectral li- 
braries and theoretical modeling are neither sufficie nt for 
verifying nor re ducing the uncertainties our models ( |Con-| 
roy et al.|2009 ) . This is largely due to the great difficulty 



in modeling post-main-sequence evolution, including the 
envelopes of TP-AGB st ars, though the UV may pro 



vide further constraints (Buzzoni & Gonzalez-Lopezlira 
20081. We are optimistic that improved characteriza- 
tion of the mid-IR colors of the TP-AGB can be incor- 
porated into SED fitting, though our calculations have 
uncertainties perhaps on the order of a factor of two due 
to uncertainties in the ensemble colors of the TP-AGB 
populations at different ages. With refinement we antici- 
pate that incorporating the mid-IR into multiwavelength 
analysis of SEDs will provide the strongest constraints on 
the star formation histories of galaxies. 



suolo et al.||2007[ ). But the model presented here implies 
that the mid-IR provides the total mass in stars formed 
in windows stretching back 1.5 Gyr in cosmic time. As 
a result, such observations must be used with care when 
constraining the star formation rate density of the uni- 
verse at z < 2, or when considering whether variations in 
the initial mass function are wa rranted by the data (e.g. 
[Dav^[2008{ [Wilkins et aLl[2008| ). 

T he detection of galaxies in trie mid-infrared over most 
of a Hubble time has helped change our view of galaxy 
assembly, and the determination of star formation rates 
associate d with that assembly has remained a difficult 
task (seelChary & Elbaz|2001 



Calzetti et al |2007| [Saliin 



et al. 2009[ |Rieke et aL||2UU9 r Perhaps the mosT lm- 
portant implication ot this Letter is that modeling the 
TP-AGB has allowed us to derive the relationship be- 
tween mid-IR luminosities and star formation rates from 
"first principles" for the first time. With such models, it 
should now be possible to more accurately constrain the 
detailed history of star formation in the universe back to 
early times. 
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